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ABSTRACT 


Conical  shaped  vehicles  with  high  ballistic  coefficients  are  subjected  during  re-entry,  to 
intense  fluctuating  air  pressures  which  result  in  high  internal  sound  and  vibration  levels. 
As  a  result  of  these  high  vibration  levels,  internal  components  can  be  adversely  affected. 
In  order  to  insure  reliability,  vibration  test  criteria  are  predicted  based  on  limited 
experimental  data  at  high  Mach  numbers  and  generally  result  in  conservatism  of  the 
design.  This  program  using  prior  experimental  data  and  high  Mach  number  data  obtained 
in  the  present  contract,  reviews  existing  analytical  expressions  necessary  to  define 
aeroacoustic  loads  associated  with  high  beta  re-entry  vehicles  and  recommends  analytical 
expressions  for  transitional,  fully  developed  turbulent,  separated  and  base  flow. 
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SECTION  I 

INTRODUCTION 


During  re-entry  into  the  earth's  atmosphere,  vehicles  are  subjected  to  intense  fluctuating 
aerodynamic  pressures .  As  a  result,  significant  vibration  response  can  occur,  affecting 
the  performance  of  components  within  the  vehicle.  In  order  to  accurately  assess  this 
problem,  .t  is  necessary  to  define  the  nature  of  the  acoustic  environment  for  various 
aerodynamic  flows  with  vehicles  during  re-entry.  This  has  become  more  important  as 
mission  requirements  have  resulted  in  higher  velocities  at  relatively  lower  altitudes. 

The  purpose  of  the  program  undertaken  in  the  present  study  was  to  measure  the  aeroacous- 
tio  environment  or.  the  surface  of  a  smooth,  conical  body  and  use  the  results  to  upgrade 
current  analytical  formulations.  The  experiments  were  conducted  at  Mach  4,  8,  and  10 
for  several  values  of  free  stream  Reynolds  number,  angle  of  attack  and  three  bluntness 
ratios . 

For  the  various  tonne'  coi—itions  and  model  attitudes^uctuating  pressure  measurements 
were  made  in  fully  turbulent  flow  and  at  the  base  of  the  teBt  vehicle.  Two  other  areas 
investigated  were  pressure  fluctuation  in  the  transition  region  where  flow  changes  from 
laminar  to  turbulent  flow  and  separated  flow  resulting  from  angles  of  attack  equal  to  the 
cone  half  angle  or  the  test  vehicle,  A  .ulficient  amount  of  acoustic  sensors  were  located 
on  the  vehicle  to  define  fluctuating  pressure  levels  acting  along  the  vehiclo. 


SECTION  II 

CONCLUSIONS  AND  RECOMMENDATIONS 


"I 


The  following  formulae  are  recommended  to  describe  the  acoustic  magnitude,  spectra 
and  cross  Spectra  and  cross  speotra  coefficients  for  various  flow  conditions: 

1.  ATTACHED  TURBULENT  BOUNDARY  LAYER 

Magnitude 
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Cross  Spectra  coefficients 
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(The  above  coefficients  are  expressions  for  the  cross  spectra  Eqs.  16,  17). 
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2.  TRANSITIONAL  BOUNDARY  LAYER  FLOW 
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4.  BASE  FLOW 
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No  significant  angle  of  attack  effects  were  observed  (pressure  spectra  or  acoustic  magnitude) 
on  pressure  fluctuation  measurements  made  in  fully  developed  turbulent  flow,  however,  the 
transition  region  on  the  test  vehicle  changed  location  as  predicted. 

The  increase  in  the  RMS  sound  pressure  level  («10  dB)  near  the  end  of  transition  evident 
in  the  data  is  similar  to  the  increase  in  the  velocity  profile  exponents,  n,  measured  by 
Martellucci*.  This  trend  suggests  that  the  velocity  profile  exponent  as  a  parameter  in  de¬ 
fining  a  single  equation  for  the  variation  of  acoustic  intensity  in  the  transition  region  as 
well  as  fully  developed  turbulent  flow. 

An  equation  has  been  developed  defining  the  peak  acoustic  magnitude  in  transional  flow 
based  on  data  obtained  for  free  stream  Mach  numbers  of  4,  8  and  10 .  This  Equation  46  is 
similar  to  that  defined  for  turbulent  flow  (3  or  6)  differing  only  in  the  coefficient  of  local 
Mach  number.  An  assumption  that  these  two  equations  are  both  applicable  at  low  Mach 
numbers  implies  convergence  of  Prms/qe  for  transitional  and  fully  turbulent  aerodynamic 
conditions  as  Me  approaches  0 .  Quantitative  effects  of  transition  at  low  Mach  numbers 
should  be  the  subject  of  further  investigations . 

The  effects  of  separated  flow  (attained  by  placing  the  model  at  an  angle  of  attack  equal  to 
its  half  cone  angle)  on  pressure  fluctuations  on  the  vehicle  were  significant.  Of  significance 
were  the  high  fluctuating  pressures  on  the  windward  side  of  the  vehicle  which  remain  high 
downstream  from  the  transition  region,  and  the  low  pressure  levels  on  the  leeward  side  of 
the  vehicle  in  the  separated  flow  region  which  result  in  a\«  15  dB  circumferential  variation 
in  fluctuating  pressures  over  a  large  portion  of  tho  vehicle.  This  will  result  in  excitation 
of  the  higher  shell  harmonics,  ’’hoop  modes",  thereby  inducing  responses  in  the  vohicle 
which  should  be  considered  in  design  analysis . 

Base  pressure  fluctuation  measurements  wore  higher  than  expected.  Appropriate  adjustments 
were  made  to  analytical  expressions  which  were  initially  verified  usiivt  subsonic  data. 


'Martellucci,  A.,  "Effects  of  Mass  Transfer  on  Hypersonic  Turbulent  Boundary  Layer 
Properties,"  A1AA  Journal ,  Voi.  10,  No.  2,  February  11)72,  pp.  1S1-1S7. 


SECTION  III 

PRIOR  EXPERIMENTAL  DATA  REVIEW 

1.  ATTACHED  TURBULENT  BOUNDARY  LAYER  FLOW 


a.  Acoustic  Magnitude 

The  overall  acoustic  magnitude  of  pressure  fluctuations  in  attached  turbulent  flow 
from  various  experimenters  is  presented  in  Figure  1.  This  summary  is  a  representative 
cross-section  of  conditions  under  which  acoustic  intensity  measurements  were  made  (flight, 
tunnel  wall,  tunnel-model).  Data  were  not  included  if  flow  conditions  were  suspect,  such  as 
Lewis  and  Banner's  X-15  data  which  probably  contains  effects  of  shocks  off  the  fin  leading 
edge,  and/or  fuselage  buffeting.  Most  experimenters  have  found  that  the  acoustic  magnitude 
(Prms)  is  best  presented  as  normalized  by  the  free  stream  dynamic  pressure  (q«). 

For  application  to  flight  vehicles  the  local  dynamic  pressure  (qe)  is  a  more  rea¬ 
sonable  choice  than  q,*,,  since  it  is  more  definitive  of  local  conditions.  Data  were  therefore 
plotted  in  Figure  1  as  a  function  of  local  properties  (Me  &  qe).  No  adjustment  was  required 
to  aviously  published  tunnel  results  since  the  local  freo-stream  properties  are  equal. 

b.  Acoustic  Power  Spectral  Density 

Figures  2  and  3d)  contain  normalized  experimental  spectra  for  subsonic  and 
supersonic  flow.  As  discussed  by  Lows  on  <2)  the  frequency  spectra  of  attached  turbulent 
boundary  layer  flow  soale  on  a  Strouhal  number  basis.  Choice  of  the  proper  typical  length 
and  velocity  is  difficult.  In  general,  tho  free  stream  velocity  (U»)  is  chosen  as  the  velocity 
parameter  and  either  boundary  layer  thickness  (<5)  or  boundary  layer  displacement  thick¬ 
ness  (  <5*)  as  the  typical  length  parameter.  In  this  study  tho  local  velocity  and  boundary 
layer  displacement  thickness  were  used  to  normalize  data  from  various  experimenters. 

A  reasonable  collapse  of  the  data  is  obtained  by  using  these  parameters  and  a 
definite  trend  In  the  data  is  apparent.  No  significant  difference  between  subsonic  and 
supersonic  flow  is  evident.  Figure  4  presents  the  same  high  Mach  number  data  normalized 
to  the  boundary  layer  thickness,  (5 ).  A  similar  spread  of  data  is  evident.  The  trend  in 
the  magnitude  of  the  spectra  to  Increase  at  Strouhal  numbers  below  3  x  10“2  is  believed  to 
be  due  to  effects  of  tunnel  noise.  At  high  Strouhal  numbers  (>1)  a  decrease  in  spectral 
values  of  approximately  20  dB/decade  is  apparent. 

c.  Cross  Correlation  Functions 

Spatial  correlation  properties  of  a  fluctuating  pressure  field  can  be  obtained  by 
•  veral  approaches.  Whited)  presents  two  forms  of  narrow  band  correlation  which  yield 
essentially  the  same  results.  Correlation  functions  may  be  obtained  by  direct  averaging 
and  normalization  of  narrow  band  components  resulting  in  data  of  the  form  shown  in 
Figures  5  and  6.  Figure  7  illustrates  narvow  band  correlation  coefficients  obtained  using 
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Figure  4.  Normalized  Power  Spectra  of  Pressure  Fluctuation  for  M„>i  Using  5 
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an  alternate  approach  which  essentially  adjusts  the  delay  times  of  the  space-time  correla¬ 
tion  functions  to  a  maxima.  Bull (8)  presents  his  findings  in  this  form  where  he  brings  out 
the  fact  that  at  high  frequencies  values  of  £ / S  *  tend  to  lie  on  a  common  curve;  however  at 
low  frequencies  there  is  significant  spread  in  the  data,  indicating  that  the  correlation  co¬ 
efficients  tend  to  be  independent  of  frequency  in  this  region.  Therefore,  expressions 
defining  cross  correlation  coefficients  require  further  adjustment  at  the  low  frequencies. 
This  can  be  accounted  for  by  utilizing  an  exponential  function  based  on  a  normalized  sep¬ 
aration  distance.  Figure  8  presents  asymptotic  values  of  narrow  band  axial  and  lateral 
pressure  coefficients  (  o>  =  0)  as  presented  by  Bull.  (8) 

The  broad  band  space-time  correlations  between  two  pressure  time  histories  arc 
useful  in  defining  behavior  of  the  flow  field.  Figures  9  and  10  (lateral  curves  represent 
correlation  peaks  at  t  =0)  are  typical  results  obtained  by  various  investigators.  Results 
show  the  typical  decrease  of  peaks  with  increasing  separation  distance,  characteristic  of 
decay  in  the  turbulent  eddy  magnitude  with  distance  from  the  reference  point.  Also  evident 
at  all  Mach  numbers  was  the  sharp  peaks  present  for  small  separation  distances.  As  the 
separation  distance  increases  these  sharp  peaks  become  broader  and  flatter.  This  effect 
also  occurs  in  the  lateral  or  cross  stream  direction.  BuliW  relates  this  behavior  of  the 
space-time  correlation  curves  to  the  presence  of  a  broad  spectrum  of  frequencies  at  small 
separations.  At  larger  separations  the  spectrum  of  correlation  producing  components  is 
narrow  and  dominated  by  low  frequency  components. 

d.  Convection  Velocity 

Bull(8)  defines  broad  band  convection  velocity  as  the  weighted  mean  convection 
velocity  of  a  pressure  source  In  the  boundary  layer  which  contributes  to  the  correlation. 
Convection  velooity  is  determined  by  dividing  the  separation  distance  between  two  sensors 
by  tho  time  betwoen  peaks  of  the  respective  space-time  correlation  functions.  Figure  1 1 
Is  a  summary  of  broad  band  convection  velocity  as  a  function  of  separation  distance  from 
various  experimenters.  Edge  velocity  and  boundary  layer  displacement  thickness  are  used 
as  normalizing  parameters,  Tito  variation  of  velocity  with  Increasing  spatial  distance  from 
0.5  U()  to  0,86  Ue  Is  attributed  to  outward  dispersion  of  turbulent  eddies  originating  in  the 
region  of  transition  botweon  the  viscous  sublayer  and  inner  pari  of  the  fully  turbulent  region 
of  the  constant  stress  layer. 

Narrow  band  convection  velocities,  which  define  the  frequency-velocity  relation¬ 
ship  associated  with  turbulent  eddies  progressing  with  the  flow,  from  various  investigations 
are  presented  in  Figure  12.  At  frequencies  below  5  kHz,  the  convection  velocity  approaches 
unity  which  indicates  that  the  low  frequency  components  of  pressure  fluctuations  are  associ¬ 
ated  with  turbulent  eddies  in  the  outer  region  of  the  boundary  layer.  High  frequency  com¬ 
ponents  of  the  pressure  fluctuations,  1x3 cause  of  their  low  convection  velocity,  are  the 
result  of  turbulent  eddies  near  the  wall. 

in  oorv'i-lng  Figures  11  and  12  it  becomes  apparent  that  the  broad-band  convection 
velocity  is  a  goou  representation  of  tile  narrow  band  convection  velocities  at  Stroulial 
numbers  less  than  two. 
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Figure  12.  Variation  of  Narrow  Band  Convection  Voloclty  With  Frequency 
for  Various  Separation  Distances 

2.  TRANSITIONAL  BOUNDARY  LAYER  FLOW 

Limited  experimental  data  describing  acoustic  intensity  and  spectral  distribution  of  trans¬ 
itional  flow  is  available.  No  data  were  found  defining  correlation  functions  or  convection 
velocities.  This  flow  regime  is  of  considerable  interest  because  vibration  data  measured 
on  high  ballistic  coefficient  re-entry  vehicles  wore  shown  to  bo  more  severe  during 
transitional  flow. 

a.  Acoustic  Magnitude 

Pate  and  Brown  measured  the  acoustic  magnitude  in  transitional  flow  on  a 
6  degree  conical  vehicle  at  Mach  3  and  4.  Tests  were  conducted  at  AEDC  (Tunnel  A 
facility)  where  the  stagnation  pressure  was  varied,  thereby  shifting  the  location  of  the 
transition  region  along  the  test  vehicle.  Using  this  approach,  pressure  fluctuations  associ¬ 
ated  with  laminar,  transition  and  fully  turbulent  flow  were  measured  using  a  single  sensor. 
Figure  13  gives  the  rms  pressure  fluctuations  for  the  frequency  range  from  0  to  25  klk. 
Similar  tests  were  also  performed  with  a  5  degree  cone-cylinder  by  Saunders  and 
Johnson^)  where  acoustic  measurements  in  transitional  flow  were  obtained.  11ie.se  re¬ 
sults  are  also  shown  in  Figure  13,  The  Increase  in  fluctuating  pressure  in  the  transition 
region  is  clearly  evident  at  both  Mach  numbers  for  all  tests. 

Figure  14  presents  the  available  acoustic  magnitude  measured  during  transition. 
Included  In  this  figure  data  point  obtained  during  re-entry  flight  of  a  high  ballistic 
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Figure  14.  Hoot  Mean  Square  Pm  a  sure  Fluctuations  Vs.  Free  St  mam 
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coefficient  re-entry  vehicle.  Telemetry  data  for  this  flight  wore  limited  to  I  kHz;  however, 
the  data  point  shown  in  Figure  14  was  obtained  by  extrapolating  the  available  data  using 
Houbolt's  distribution U4)  for  fully  turbulent  flow.  The  normalized  rms  pressure  out  to 

4  kHz  had  a  value  of  Pnns/q  *  0.0001S. 

b.  Acoustic  Power  Spectral  Density 

The  spectral  distribution  in  transitional  flow  are  shown  in  Figure  15  for  measure¬ 
ments  made  in  wind  tunnel  experiments.  Results  show  a  considerable  amount  of  energy  in 
the  low  frequency  range.  Figure  16  is  a  comparison  of  spectra  obtained  during  transition 
and  fully  turbulent  (Mg  >10)  flow  from  flight  tost  data(lS).  The  characteristic  high  levels 
at  low  frequencies  arc  also  evident  in  the  spectra,  although  not  os  severe  as  was  shown  in 
wind  tunnel  experiments.  Though  no  firm  conclusions  can  be  drawn  because  of  the  limited 
amount  of  data,  it  appears  that  a  considerable  amount  of  the  low  frequency  spectral  data 
may  be  attributed  to  tunnel  noise  and  electrical  effects.  Shown  In  Figure  17  arc  typical 
acceleration  spectral  densities  measured  on  a  internally  mounted  component  during  a 
re-entry  vehicle  flight  (Me  >10),  Comparing  measurements  made  during  transitional 
and  fully  turbulent  flow,  the  transition  period  is  seen  to  be  more  severe  by  approximately 

5  dll  at  peak  and  6  dB  throughout  the  remainder  of  the  spectra.  This  is  lower  than  observa¬ 
tions  made  on  external  environments  where  fully  developed  turbulent  and  transitional  flow 
differ  by  about  io  dB. 
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Figure  15.  Transitional  Flow  Spectral  Distribution  Obtained  in  Wind 

Tunnel  Experiments 


EHEOUENCY.  M/ 

Figure  16.  Comparison  of  Transition  !»nd  Turbulent  Flow  Pressure  Spectra 

Measured  In  Flight 


FREQUENCY.  H/ 

Figure  17.  Comparison  of  Acceleration  Speotra  Measured  in  Transitional 

and  Fully  Turbulent  Flow 


3.  SEPARATED  FLOW 
a.  Acoustic  Magnitude 

A  review  of  separated  flow  data  given  by  Robertson  in  Reference  1  divides  the 
separated  flow  into  various  categories,  several  of  which  are: 

(1)  expansion  induced  separated  flow 

(2)  separated  flow  upstream  of  compression  corner 

(3)  flare  induced  separation 

(4)  shock  wave  interaction 

The  following  is  a  brief  summary  of  data  presented  In  the  above  reference. 

”tgure  18  represents  the  overall  acoustic  magnitude  normalized  by  the  free  stream  dynamic 
'.treasure  as  a  function  of  local  Mach  number  for  various  separated  flow  environments  down** 
.stream  from  expansion  corners.  Such  expansion  Induced  separated  flow  is  representative 
of  regions  aft  of  eonc-cyltnder  junctions,  rearward  facing  steps  and  in  the  near  wake  of 
boat  tail  configurations.  Tolerance  brackets  shown  In  the  data  represent  variations  due  to 
nonhomogeneous  flow  within  the  region  of  constant  static  pressure  rather  than  scatter  in  the 
measurements. 
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Figure  18.  Overall  Acoustic  Magnitude  Vs.  Local  Mach  Number  for  Expansion 

Induoed  Separated  FlowU) 

Data  representative  of  flare  Induoed  separation  and  separated  flow  upstream  of 
compression  oornors  is  given  in  Figure  19,  No  analytical  expressions  have  been  developed 
for  this  type  of  flow  since  insufficient  data  in  regards  to  local  flow  conditions  is  available. 
In  general  the  overall  pressure  levol  appears  to  Increase  slightly  with  increasing  Mach 
number.  This  may  be  due  to  the  Improper  choice  of  ordinate  variable. 

The  oscillation  of  a  shock  wave  produces  intense  fluctuating  pressures  for  the 
regton  in  dose  proximity  to  tire  shook  wave.  Included  In  Figure  20  as  examples  of  shock- 
wave  oscillation  mid  reattachment  data  are  measurements  by  Chyu  and  iianly.  (*) 

Speaker  and  Allman(^)  have  also  defined  overall  intensity  levels  for  the  separated  flow 
reattaohment  as  well  as  shock  Interaction.  These  also  appear  in  Figure  20, 

b.  Acoustic  Power  Spectral  Density 

A  summary  of  available  data  for  separated  flow  power  spectra  obtained  in  the 
homogeneous  region  of  compression  is  given  in  Figures  21  and  22*  Heohtoin'^  performed 
a  study  (at  Mach  2)  where  various  cliarueterlstlo  lengths  were  used  in  normalizing  the 
frequency  parameter  along  with  the  mean  local  velocity.  Characteristic  lengths  such  as 
calculated  displacement  thickness  in  tile  separated  flow  region,  calculated  displacement 
thickness  associated  with  the  zero  velocity  streamline,  and  distances  of  the  measurement 
behind  the  shock  wave  were  used.  Of  these  parameters,  the  distance  of  the  measurement 
behind  the  shock  was  the  most  appropriate  length  parameter  for  frequency  scaling.  This 
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Figure  19.  Overall  Acoustic  Magnitude  Vs.  Free  Stream  Mach  Number  for 
Separated  Flow  Upstream  of  Compression  Comers U) 

length  parameter  was  also  shown  to  be  most  appropriate  in  the  normalization  of  the 
amplitude  scales  of  the  power  spectral  densities.  Other  parameters  used  to  normalize  the 
spectra  amplitude  in  conjunction  with  the  length  parameter  wero  local  static  pressure, 
free  stream  dynamic  pressure  and  the  square  of  the  Mach  number.  Of  these,  Urn  free 
stream  dynamic  pressure  was  found  to  collapse  the  data  best.  Figure  23  compares  this 
normalization  technique  to  tliat  used  by  Cliyu  and  Hanly,  who  employed  the  boundary  layer 
displacement  technique  in  the  attached  flow  region  prior  to  the  shock  front,  free  stream 
dynamic  pressure  and  velocity  as  normalization  parameters.  Both  normalizing  techniques 
appear  to  collapse  the  data  well  with  slightly  better  collapse  of  data  using  attached 
boundary  layer  displacement  thickness.  Figure  24  is  a  comparison  of  spectral  data  for 
expansion  induced  separated  flow  and  compression  induced  flow  obtained  by  Speaker  and 
Ailmanw)  at  Mw  -  3,40.  Though  the  curves  arc  several  orders  of  magnitude  apart,  their 
shapes  are  similar.  The  difference  in  spectral  level  is  probably  due  to  normalizing 
parameters. 

Robertson (1)  presents  spectral  characteristics  of  shock-wave  oscillation  for  both 
two  and  three-dimensional  protuberances  (given  in  Figure  25),  The  spectrum  shows  a  sig¬ 
nificant  amount  of  energy  in  the  low  frequency  range  and  a  steep  roll-off  of  approximately 
8  dB/octave  for  a  Strouhal  frequency  range  of  1  x  10“2  to  2  x  10-1.  Speaker  and  Allman 
have  also  conducted  experiments  examining  the  spectral  distribution  in  the  vicinity  of  a 
two-dimensional  shock  wave  Impinging  upon  the  turbulent  boundary  layer,  as  well  as  in  the 
region  of  flow  rcattachmcnt,  with  similar  results. 
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figure  20.  Acoustic  Magnitude  Distribution  Representing  Shock 
Sejiaration  aad  Roattachment 


1*  Normalized  (Using  b)  Pressure  Spectra  Representative  of  Separated  Flow 
(Prel.  Data  Taken  from  Ref.  16  With  Permission  From  NASA  Ames) 


Figure  22.  Normalized  (Using  S-Xs)  Pressure  Spectra  Representative  of  Separated  Slow 
(Preliminary  Data  Taken  From  Ref.  16  With  Permission  From  NASA  Ames) 


Comparison  of  Normalization  Technique  for  Separated  Flow  Spectra 


COMPRESSION 


EXPANSION  INDUCED 


c.  Cross  Correlation  Functions 


Normalized  narrow  band  longitudinal  cross  spectra  for  separated  flow  are 
presented  in  Figure  26  for  the  Mach  No.  range  of  1.6  to  2.5(4),  These  cross  spectra  are 
similar  to  those  obtained  for  fully  turbulent  flow  and  are  represented  as  damped  sinusoids. 
The  rate  of  decay  is  approximately  50%  greater  in  separated  flow  when  compared  to  full 
turbulent  flow  and  a  slight  increase  in  the  rate  of  decay  was  observed  as  the  Mach  number 
was  decreased.  These  effects  become  more  pronounced  in  the  coherence  functions  shown 
in  Figure  27. 

Broad  band  space  correlation  functions  for  the  above  experiment  are  presented  in 
Figure  28.  Similar  to  that  observed  from  narrow  band  analysis  was  the  fact  that  correla¬ 
tion  in  the  attached  flow  was  significantly  higher  than  that  in  the  separated  flow. 

Shown  in  Figure  29 U)  is  the  degree  of  correlation  which  exists  between  the  region 
under  a  detached  shock  wave  and  adjacent  regions  of  attached  and  separated  flow.  Weak 
correlation  is  noted  for  adjacent  regions  with  strong  correlation  existing  within  the  shock 
region. 


d.  Convection  Velocity 

For  separated  flow,  the  narrow  band  convection  velocity  varies  from  0,2  U0 
(at  low  wave  numbers)  to  approximately  unity  at  high  wave  numbers,  thus  convective 
velocities  derived  from  broad  band  analyses  (cross  correlation  functions)  aro  not  repre¬ 
sentative  of  all  frequency  components  making  up  the  pressure  field  as  was  the  case  for 
attached  turbulent  flow.  Figure  30(4)  presents  typical  separated  flow  (compression  corner 
induoed  separation)  narrow  band  convection  velocities  for  various  separation  distances  at 
Mach  2.  Figure  30  compares  the  narrow  bum)  convection  velocities  for  separated  and 
attached  flow.  IX) tailed  studies  of  the  narrow  band  convection  velocity  In  separated  flow 
were  conducted  by  Hechtlen.  U6)  Typical  results  of  his  studies  are  presented  in  Figure  31 
where  the  frequency  scale  (normalized  by  distance  from  separation  point  and  local 
velocity,  Figure  32  unnormalized)  has  been  divided  into  three  regions  in  order  to  describe 
tlie  mechanisms  attributed  to  the  pressure  field  in  separated  flow.  Low  frequency  com¬ 
ponents  of  tiie  pressure  field  below  Strouhai  numbers  of  0,15  are  a  result  of  slow  re¬ 
circulating  fluid  in  the  region  underlying  the  high  speed  separated  shear  flow.  (Significant 
scatter  in  data  is  attributed  to  small  phase  angles  used  to  calculate  convection  velocities  at 
low  Strouhai  numbers.)  For  Strouhai  numbers  between  0,15  and  1.1  measured  fluctuations 
are  due  to  convoeted  patterns  in  the  outer  regions  of  the  separated  boundary  layer.  Above 
a  Strouhai  number  of  1,1,  pressure  fluctuations  art*  due  to  eddy  Mach  wave  radiation 
originating  in  a  limited  region  immediately  behind  the  flow  separation  shock  wave. 

4.  BASE  FLOW 

As  was  the  case  for  transitional  flow,  only  limited  data  were  available  defining  the  pressure 
sjjoetra  and  overall  magnitude  for  base  flow.  No  data  defining  correlation  functions  or  con¬ 
vection  velocities  were  found. 


Figure  27.  Coherence  Function  lor  Attached  and  Separated  Flow (4) 
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Figure  29.  Correlation  of  Pressure  Fluctuations  Between  Shock  Wave  ;uul 
Adjacent  Attached  and  Separatod(l)  Flow 

a.  Base  Flow  Acoustic  Magnitude  and  Spectra 

Eldrecl(18)  published  experimental  results  at  subsonic  flow  (V  ■  ti8  to  332  i'psi  «»f 
random  pressure  fluctuations  for  base  flow.  Contained  in  die  table  insert  of  Figure  33  arc 
normalized  overall  pressure  levels  for  two  locations  on  a  circular  base  at  several  Mach 
numbers.  The  apparent  trend  for  the  normalized  overall  intensity  is  to  increase  as  one 
moves  away  from  the  center  of  the  plate. 

Also  shown  In  Figure  93  is  the  typical  spectral  distribution  measured  during  these 
tests.  Data  was  scaled  to  a  Stroulml  number  using  the  base  diameter  as  the  characteristic 
length  and  free  stream  velocity  as  the  associated  velocity.  As  expected,  a  significant 
amount  of  energy  is  associated  in  the  low  frequency  region. 
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1.  Narrow  Band  Convection  Velocities  Using  Separation  Dis Lance  as  Normal 
Preliminary  Data  Taken  from  Ref.  IG  With  Permission  from  NASA  Ames) 


2000 


trc  32.  Typical  Uanormalized  Narrow  Band  Convection  Velocity  for  Separated  Flow 
(P ret* mittary  Data  Taken  from  Ref.  IS  With  Permission  From  NASA  Ames) 
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Figure  33.  Typical  Spectral  Distribution  in  Base  Flow  (Subs^io 


SECTION  IV 


PRIOR  PREDICTION  METHODS  REVIEW 

This  section  briefly  discusses  existing  mathematical  models  of  aero-acoustic  environments 
for  various  type  flows  expected  during  the  re-entry  phase  of  high  ballistic  coefficient  ve¬ 
hicles  .  Extensive  analytical  and  experimental  studies  have  been  performed  by  several  in¬ 
vestigators  for  defining  acoustic  properties  during  turbulent  flow.  Limited  prediction 
techniques  have  been  published  for  separated  and  base  flow.  No  mathematical  models  for 
transitional  flow  were  found  to  exist  in  the  literature. 

1 .  ATTACHED  TURBULENT  BOUNDARY  LAYER  FLOW 


a.  Acoustic  Magnitude 

It  is  not  the  intent  of  this  study  to  summarize  the  historical  background  of  pres¬ 
sure  fluctuations  in  turbulent  flow,  since  that  subject  has  been  the  source  of  theoretical 
and  experimental  investigations  for  many  years .  What  will  be  reviewed  are  present  me¬ 
thods  used  in  defining  aero-acoustic  environments  during  turbulent  flow.  Houbolt  (14)  per 
formed  a  study  assuming  that  eddy  velocity  is  proportional  to  free-3tream  velocity  and  the 
local  mean  density  in  the  region  of  maximum  noise  generation  was  the  significant  variable 
governing  noise  production.  Using  Euler's  equation  for  flow  and  Crocco's  equation  for  the 
temperature  velocity  relation,  the  rms  pressure  levels  (Equations  l  &  2)  as  a  function  of 
Mach  number  and  dynamic  pressure  (or  free  stream  pressure  P)  is  defined  as: 


P 

rms 


0.007 

- k  q 

1  +0.012  M 


or 


„  _  0.0049  M2 

Prms  “  2  P 

1  +0.012  M 


(1) 


(2) 


Houbolt  in  unpublished  work  has  further  refined  his  above  expression  to  include  temperature, 
gas  density  and  a  recovery  factor  to  obtain: 


(2) 
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(4) 


VTe  =  1+ri 
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(5) 


Similar  studies  performed  by  Lowson(2)  assuming  an  adiabatic  wall,  constant  static  pres¬ 
sure  through  the  boundary  layer,  and  Crocco's  relation  resulted  in  the  following  equations 
for  the  overall  fluctuating  pressure  intensity: 


0.006 


rms 


1  +  0.14  M 


(6) 


0.0042  M2 

prms  “  - -  P  (7) 

1  +  0.14  NT 

Ailman(ia)  developed  an  expression  for  the  pressure  intensity  based  on  a  review 
of  test  results  for  the  Mach  number  range  from  0,6  to  5.  This  expression,  given  in 
Equation  8,  must  be  used  with  some  reservations  at  Mach  numbers  greater  than  five. 


Figure  34  gives  the  pressure  intensity  normalized  by  dynamic  pressure  as  a  func¬ 
tion  of  Mach  number  for  the  various  investigators  diseussod.  Of  the  various  prediction 
techniques,  Lowson's  equation,  because  of  its  simplicity  and  good  agreement  with  experi¬ 
mental  data  was  found  to  be  most  representative. 

prms  =  x  10"4  V6  +  (M«  "  4)2] 


b.  Acoustic  Power  Spectral  Density 
(21) 

Bull  and  Willis*'  semi-empirical  approach  defined  the  spectral  relation  in 
turbulent  flow  as: 
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‘e  34.  Comparison  of  Overall  Acoustic  Magnitude  Prediction 
Methods  to  Experimental  Data 


Later,  Bull (^2)  using  miniature  pressure -sensing  elements  in  subsonic  wind 
tunnel  tests  and  auto -correlation  functions  developed  the  following: 
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where 


«0  =  W^VUoo 


(14  20) 

Using  experimental  data  and  engineering  logic,  Houbolt'  ’  derived  a  spectral 
relation  for  turbulent  flow  as: 


25  2 
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rms 
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Lowson^  using  subsonic  and  supersonic  data  developed  an  empirical  Equation 
(12)  which  introduced  the  boundary  layer  thickness  as  a  representative  length  parameter. 
This  is  given  by  the  following: 


0(w)  = 


prms 
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Robertson^  using  recent  measurements  at  supersonic  speeds  developed  Equation 
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where 


u;Q  *  o.5  U«/&‘ 
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This  Equation  13  does  give  the  best  fit  to  experimental  data  as  shown  in  Figure  35; 
however,  Houbolt's  Equation  11  is  recommended  for  general  use  because  of  its  simplicity. 
Since  no  one  characteristic  length  (s  or  5*)  gives  better  collapse  of  the  data,  Lowson's 
Equation  is  also  recommended  for  use  if  boundary  layer  thickness  is  known.  Figure  36 
compares  Lowson’s  Equation  with  typical  spread  in  test  data. 

c.  Cross -Correlation  Functions 


Cross-correlation  and  cross-spectrum  functions  pertinent  to  the  aero-acoustic 
environment  are  related  through  their  Fourier  transforms  as  follows: 


<t>  ( £»  n>  u>) 


i  art 

R  ( £,  v ,  t)  e 


dt 
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and 
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dor 
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The  cross -spectrum  can  also  be  expressed  in  terms  of  the  co-spectral  density 
and  quad  spectral  donsity; 

*  (*.»?»«)  =  0r  ( I.  >>.  «)  -  j  0  Ui9.«)  (15) 
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whore  the  co-spectral  density  is  used  in  defining  characteristics  of  tho  fluctuating  pressure 
field.  Lowson(2)  and  RobertsonU)  using  data  obtained  by  Bull  (8)  defined  the  cross -cor  rela¬ 
tion  function  for  turbulent  flow  assuming  that  the  cross-power  spectral  density  is  a  function 
only  of  the  separation  distances.  These  functions  are  of  tho  form: 
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From  Bull’s  data  the  following  expressions  for  the  correlation  coefficients  are  given  as: 
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Figure  35.  Comparison  of  Normalized  1*0 wo r  Spectra  Prediction  Methods 

to  Experimental  Data 


Figure  36.  Comparison  of  Lowson's  Normalized  Power  Spectra 
(Using  6  )  to  Experimental  Data 


These  correlation  coefficients  may  also  be  expressed  in  terms  of  boundary  layer  displace 
ment  thickness  as: 


A,,(n,  cj  )  =  exp 
(Figure  36) 
A|($,  <*>)  ®  exp 

(Figure  37) 


(Figure  39) 
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(Figure  3D) 


(20) 


(21) 


Lewson^  suggests  that  when  using  the  separable  form  of  the  correlation  coefficient  in 
response  analysis,  the  correlation  area  is  underestimated  by  n  factor  of  */2  and  should  be 
adjusted  accordingly.  This  factor  stems  from  the  fact  that  the  coefficients  will  remain 
constant  along  straight  lines  on  the  surfaces,  forming  a  diamond  pattern  surrounding  the 
origin  which  physically  is  unreasonable.  An  elliptical  form  would  be  more  realistic. 
Therefore  a  ratio  of  areas  is  required,  leading  to  a  ir/2  factor.  Correlation  functions 
defined  by  Crocker^23)  and  White  are  of  similar  forms  to  those  previously  discussed 
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Figure  37 .  Comparison  of  Frequency  Dependent  Narrow  Bancs  Dong  apace 
Correlation  Coefficient  with  Analytical  Expression 
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Figure  38.  Comparison  of  Frequency  Dependent  Narrow  Band  Lat.  Space 
Correlation  Coefficient  With  Analytical  Expression 


except  the  coefficient  of  0.72  in  Equation  18  is  specified  as  2.0.  Subsequent  published 
data  suggests  that  a  coefficient  of  0.72  is  more  reasonable  (Figure  38).  Maestrello(^) 
suggests  the  cross-correlation  coefficient  to  be: 
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(22) 


Using  Maestrello's  data  at  Me  =  0.6  results  in  the  following  expression  for 

-1.71*1  -3 . 9 1 7 1  -iwi/u 

6  7"e  T" e  c 

which  suggests  a  significantly  higher  decay  rate  than  that  defined  by  Bull  (Equations  18,  19). 

Houbolt  (Equation  23)  defined  a  similar  general  expression  for  the  cross-correla¬ 
tion  function  as  Maestrello,  however  specific  longitudinal  exponential  values  were  not 
presented. 


A  ( S, »? »  <j) 
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where  C  “  1,  a  -  to  be  defined,  varies  with  Mach  No. 


(23) 


Neither  Houbolt  nor  Maestrello  defined  exponent  docay  torms  as  functions  of 
tu/Uc.  This  approach  appears  acceptable  in  the  practical  design  of  high  beta  vehicles, 
duo  to  the  low  frequency  range  of  Interest  (f  L  2khz)  and  high  convective  velocities,  caus¬ 
ing  the  exponential  term  (e’H^c)  to  approach  unity. 

d.  Convection  Velocity 

(2) 

Using  experimental  data,  Lowson  defines  the  convection  velocity  normalized  to 
tlio  free  stream  velocity  as; 
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where  the  broad  band  convection  velocity  corresponds  to  ~  * 
resulting  in 
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A  corresponding  equation  for  the  broad  band  convection  velocity  as  a  function  of  boundary 
layer  displacement  thickness  would  be 

=0.8- 0.5  exp  (-0.4*/S#)  (26) 

U.O 

2.  SEPARATED  FLOW 

A  comprehensive  treatment  of  separated  flow  and  oscillating  shocks  in  terms  of  pre¬ 
diction  methods  is  given  by  Robertson(l) .  This  section  briefly  summarizes  those  results. 

a.  Acoustic  Magnitude 


For  expansion  induced  separated  flow,  Equation  27  is  recommended: 
0.045 

P  =  -  Q 

rms  2  00 

1+Me 


(27) 


The  expression  given  by  Equation  27  is  compared  with  experimental  data  in  Figure  40 . 

This  expression  is  similar  in  form  with  that  proposed  for  attached  turbulent  flow.  At  this 
time  no  general  expression  is  available  for  defining  fluctuating  pressure  lovels  in  the 
region  upstream  of  compression  comers,  since  insufficient  experimental  data  is  available. 
The  data  as  shown  in  Figure  19  appear  to  increase  with  free  stream  Mach  number.  It 
would  seem  reasonable  that  local  flow  conditions  would  play  an  important  part  in  defining 
an  expression  for  the  overall  fluctuating  pressure. 

At  the  point  of  flow  separation,  the  flow  intermittently  fluctuates  between  a 
separated  flow  and  attached  flow  condition.  This  represents  an  alternating  unbalance 
between  the  large  pressure  rise  through  the  shock  wave,  exceeding  that  required  for 
separated  flow.  Hence,  extremely  large  fluctuating  pressures  can  result  from  this  con¬ 
dition.  Data  is  limited  for  this  environment  and  no  empirical  equation  has  been  uncovered 
in  the  literature.  Using  limited  data  of  Referenced  an  expression  for  fluctuating  pressures 
at  the  point  of  separation  is  proposed  as: 

Prms  a  — -  (28) 

1  +0.5  MJ* 

M) 

Figure  41  compares  available  data  with  Equation  28.  Chyu  and  Manly  present 
fluctuating  pressure  data  normalized  by  the  free  stream  dynamic  pressure  at  the  point  of 
shock  reattachment.  These  levels  are  significantly  higher  than  at  Uie  point  of  separation. 
However,  insufficient  data  is  available  to  derive  a  representative  expression  for  fluctuating 
pressure  in  the  region  of  shock  wave  rcattachmcnt. 
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Figure  *0-  Comparison  of  Overall  Acoustic  Intensity  for  Separated  Flow  With  Analytical  Predictions 
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Figure  41 .  Comparison  of  Overall  Acoustic  Mjignitude  et  the  point  of  flow 
Separation  With  Analytical  Predictions 

b.  Acoustic  Power  Spectral  Density 

Using  data  obtained  by  Hechtien^1^  for  Urn  homogeneous  region  of  compression 
corners,  Robertson  defined  the  normalized  power  spectra  as: 
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Equation  29  apjjoars  to  be  a  reasonable  expression  for  expan si  n  induced  separated  flew 
even  though  it  is  derived  based  on  datn  taken  in  compression  corners,  as  indicated  in 
Figure  24.  Figure  42  compares  this  expression  with  typical  test  data . 

Equikloc  32  defines  the  power  spectral  density  distribution  for  shock  wave  oscillation  in 
homogeneous  flow.  Hobortson  derived  tills  expression  by  combining  the  spectral  distribu¬ 
tions  for  separated  flow  (Equation  29)  and  shock  wave  oscillation,  Equation  30,  (in  the  ab¬ 
sence  of  viscous  flow)  with  an  expert  mentally  derived  weighting  function  (Kj).  Flow  param¬ 
eter  used  in  these  expressions  are  for  conditions  upstream  of  the  shock  wave. 


S3 


Separated  Flow  Powe 


—  Overall  level  of  shock  oscillation  peak  corresponding  to 
•he  mean  location  of  the  shock  wave. 


Overall  level  of  homogeneous  separated  flow. 


Subscripts :  SW  -  Shock  wave 

S  -  Sc |«ira tod  flow 

Superscripts:  I  -  Absence  of  viscosity 
H  -  Homogeneous  flow 

= 1 x 10"2 
SW 

A  comparison  of  the  predicted  spectral  distribution  for  shock  wave  oscillation  is  presented 
in  Figure  43. 

o.  Cross-Correlation  Function 

Chyu  and  HunlyH)  data  indicate  tint  decay  of  the  normalized  cross  spectra  for 
separated  flow  is  exponential  at  high  values  of  w  t  /Ue,  similar  to  that  observed  for 
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Figure  43.  Comparison  of  Shock-Wave  Oscillation  Spectral  Distribution 

to  Analytical  Prediction 

attached  flow.  Thus  as  a  first  approximation,  the  normalized  longitudinal  co-spectra 
may  be  ropi*esented  by: 
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However,  the  associated  decay  rates  were  found  to  vary  for  each  tost  Mach  number 
(1.6  to  2.5),  suggesting  Unit  the  decay  rate  may  be  dependent  on  Mach  number  as  well 
as  local  flow  conditions.  This  suggests  that  the  above  functions  (Equations  33,  34) 
should  be  expressed  In  the  following  form: 
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No  data  or  expressions  defining  the  transverse  spatial  characteristics  of  shock  induced 
fluctuating  pressures  were  found  in  tlie  literature. 

3.  BASK  FLOW 

Houbolt's  technique^®)  was  the  only  method  available  for  predicting  base  pressure 
fluctuating  characteristics.  This  method  is  based  upon  knowledge  of  tlie  characteristic  of 
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aerodynamic  turbulent  flow  and  the  application  of  engineering  logic  rather  than  experimental 
data.  Houbolt's  formulation  is  summarized  in.  the  following  sections. 


a.  Acoustic  Intensity 

The  pressure  fluctuations  in  base  flow  are  a  function  of  the  strong  shear  flow  along 
the  wake  cone  and  a  complex  recirculation  flow.  Houbolt  assumed  that  the  mean  density  at 
some  effective  location,  for  example  point  A  in  Figure  44,  has  governing  control  on  base 
fluctuating  pressures.  Assuming  that  the  eddy  velocities  are  proportional  to  the  mean  shear 
flow  velocity,  a  simple  form  for  the  rms  pressure  level  is  expressed  in  terms  of  density 
and  pressure 


resulting  in  a  general  expression  for  overall  sound  pressure  level. 
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Figure  44.  Typical  Base  How 


58 


Tentative  values  have  been  established  for  c  and  r  with  c  chosen  to  agree  with  data 
obtained  in  Reference  2.  A  value  of  0.9  was  chosen  for  the  recovery  factor  assuming  that 
a  strong  shear  and  minimal  cooling  occurs.  The  resulting  expression  for  base  fluctuations 
is  therefore: 


1  +  0.18  M, 


2  Pb 


b.  Acoustic  Power  Spectral  Density 


In  order  to  arrive  at  the  base  spectral  density  at  a  given  point  a  correlation  function 
of  the  form 


2  L 

R(t)  =  P  e 

rms, 
b 


is  assumed,  whore  V-j-  is  an  effective  transport  velocity  and  L  is  the  scale  of  turbuleneo 
defined  as: 


[”  B(x) 

J  p2 

*'o  rms, 


dx.  (x  =  VT  t) 


The  spectral  function  consistent  with  tiro  above  defined  correlation  function  is: 
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rms. 
b 


l+l  —  ) 
v  U 

o 


Assuming  that  the  scale  of  turbulence  is  on  tlie  order  of  the  base  radius  and  the  transiwrt 
velocity  Is  ix) luted  to  the  transverse  component  of  tlie  flow  velocity  along  the  wake  boundary 
tlie  following  spectral  function  results 
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where  S  is  the  apparent  wake  cone  surface  length.  Figure  45  compares  the  above  axprea- 
sion  with  experimental  data  obtained  from  Reference  7. 

c.  Cross  Correlation 


No  experimental  data  or  strict  derivation  of  the  correlation  function  was  found. 
Houbolt  has  assumed  the  following  function: 


R(t)  =  P. 


rms, 


(43) 


VT  =  Ub  sin  9 


where 


SECTION  V 


SUMMARY  OF  TEST  RESULTS 


This  section  represents  a  condensation  of  Section  n  (Volume  II)  of  this  study.  Fluctuating 
pressure  data  were  obtained  for  transitional,  turbulent,  separated  and  base  flow  at  various 
angles  of  attack  ( a),  Mach  numbers  and  nose  radii.  Several  acoustic  sensors  were  located 
on  the  surface  of  the  forward  section  of  the  test  model  (0  to  8  inches)  in  laminar  flow  to 
measure  tunnel  tare  noise.  This  term  is  used  since  in  addition  to  tunnel  noise,  some 
radiated  noise  contribution  from  turbulent  flow  downstream  is  also  measured  by  these 
sensors.  For  various  tunnel  conditions  especially  at  the  higher  Mach  numbers,  these 
sensors  measured  significant  fluctuating  pressures  which  tended  to  mask  measurements 
made  in  fully  developed  turbulent  flow.  However,  it  was  possible  to  correct  the  data  in 
order  to  obtain  realistic  approximations  to  the  true  values.  Tunnel  conditions,  model 
attitude  and  configurations  for  which  data  was  obtained  are  given  in  Table  I.  Tables  II  and 
III  define  the  free  stream  and  local  aerodynamic  properties  while  Table  IV  presents  a  sum¬ 
mary  of  acoustic  sensors  installed  in  the  test  vehicle  along  with  location  definition. 

1.  ACOUSTIC  MAGNITUDE 

Fluctuating  pressures  (from  0-20  kHz)  during  turbulent  flow  wore  obtained  from  tire  aft 
array  of  acoustic  sensors  for  tire  three  test  Mach  numbers:  4,  8  and  10.  These  pressures, 
normalized  by  the  local  dynamic  pressure  are  listed  in  Table  V,  It  was  observed  that  as 
Reynolds  numbers  were  increased,  tunnel  tare  noise  also  increased.  This  is  attributed  to 
the  inclusion  of  an  additional  stage  of  compression  necessary  for  desired  tunnel  conditions 
and  the  forward  advancement  of  the  region  of  transition  along  tho  tost  vehicle.  Actual  and 
corrected  pressure  data  at  all  Mach  numbers  are  given  in  Table  V,  which  illustrates  tho 
effect  of  tunnel  taro  noise.  Because  of  tire  uncertainty  in  tunnel  tare  noise,  it  is  recom¬ 
mended  that  the  levels  given  in  Table  V  be  considered  as  an  upper  and  lower  bound  of  tho 
true  pressure  data. 

Table  VI  summarizes  the  peak  normalized  fluctuating  pressures  under  transition  of  flow. 
The  spatial  distribution  of  the  acoustic  environment  for  transition  is  a  gradual  build-up  in 
fluctuating  pressure  reaching  a  maximum  value  toward  tho  end  of  the  transition  zone. 

Fully  turbulent  flow  is  achieved  at  distances  downstream  ranging  from  10  to  100  boundary 
layer  thicknesses  from  ousot.  Typical  fluctuating  pressure  distributions  along  tiro  vehicle 
are  shown  in  Figure  46. 

Several  tests  were  conducted  at  Mm  “8.0  with  the  test  model  at  an  attitude  of  ;j.7.2  degree 
(vehicle  half-cone  angle)  which  exposed  the  main  array  of  acoustic  sensors  to  both  wind- 
ware!  and  leeward  flow  resulting  in  a  Complete  fluctuating  pressure  map  of  the  vehicle. 
Figure  47  represents  typical  fluctuating  pressure  distributions  measured  at  these  high 
angles  of  attack.  Of  significance  In  those  runs  are  the  unexpectedly  high  fluctuating  pres¬ 
sures  on  tlie  windward  side  of  the  vehicle  which  remain  high  downstream  from  the  trans¬ 
ition  region.  Co-cxlsteut  with  these  high  levels  are  the  surprisingly  low  pressure  levels 


62 


TABLE  L,  SUMMARY  OF  NOMINAL  TEST  CONDITIONS 


‘Tunnel  doors  oj>en,  remaining  tosita  in  Tunnel  U  .vith  doors  dosed 


TABLE  II.  SUMMARY  OF  FREE  STREAM  CONDITIONS 


Facility 

^ccft 
(x  10‘6) 

Po 

(psia) 

T0 

(°R) 

u« 

(fps) 

C[oo 

(fsia) 

Poo 

(psia) 

Tqo 

(°R) 

(Lbm/ft3) 

A 

4 

1.4 

18.2 

605 

2352 

1.348 

0.120 

144 

1.28  x  10"6 

2.2 

28.4 

602 

2332 

2.109 

0.187 

143 

2.0  x  10-6 

2.8 

39.9 

597 

2332 

2.651 

0.236 

142 

2.56  x  lO-6 

3.0 

38.3 

603 

2332 

2.832 

0.252 

143 

2.70  x  10-6 

B 

8 

1.4 

286.2 

1257 

3742 

1.351 

0.0305 

92.2 

8.93  x  10-4 

2.2 

473.8 

1297 

3800 

2.212 

0.0497 

94.7 

1.41  x  10-3 

2.8 

620.9 

1321 

3836 

2.866 

0.0641 

95.9 

1 . 80  x  lO"3 

3.25 

728.2 

1320 

3835 

3.342 

0.0745 

95.7 

2.10  x  10~3 

C 

10 

1.4 

1081 

1910 

4793 

1.717 

0.0245 

95.7 

6.92  x  10~4 

2.2 

1811 

1911 

4804 

_ 

2.755 

0.0383 

93.6 

1.105  x  10’3 

on  the  leeward  side  of  the  vehicle  measured  in  the  separated  flow  region.  The  composite 
distribution  results  in  a  IS  to  20  dB  circumferential  variation  in  fluctuating  pressures 
over  a  significant  portion  of  the  vehicle  surface. 

Fluctuating  pressure  measurements  were  also  made  on  the  base  of  the  test  vehicle. 
Figure  48  presents  typical  distributions  for  three  Reynolds  numbers  at  Mach  4  and  10. 

As  Reynolds  number  Increased,  fluctuating  pressures  were  also  observed  to  increase 
proportionally. 

Hie  effects  of  bluntness  were  also  investigated  and  found  to  have  only  negligible  effect 
except  for  moving  the  transition  location  ou  the  vehicle  surface, 

2.  POWER  SPECTRAL  DENSITY 


Acoustic  spectra  for  turbulent  flow  wore  obtained  at  Mach  4,  8  and  10.  As  previously 
noted,  tunnel  biro  noise  was  appreciable  In  relation  to  turbulent  flow  data.  As  Reynolds 
number  or  Mach  number  arc  increased  tunnel  tare  also  increases.  Using  data  from 
Re*,  a  1.4  x  106  affords  the  maximum  separation  of  noise  and  true  data.  Figure  49  gives 
spectra  as  obtained  from  measured  data  with  no  consideration  to  tunnel  tare  noise. 


14.3 


Cane  Surface 


TABLE  V.  ACOUSTIC  MAGNITUDE  FOR  TURBULENT  FLOW 


Free  Stream 
Mach  No, 

Local 

Mach  No. 

prms/<le  (0-20  kHz) 

Actual  Data 

Corrected  Data 

4 

3.7 

1.62 x 10~3 

1.31  x 10"3 

8 

6.8 

1.1  x I0“3 

0.53  x  ID'3 

10 

8.1 

1.27  x  10~3 

. . . . . . 

0.67  x  10-3 

TABLE  Vi.  PEAK  ACOUSTIC  MAGNITUDE  IN  TRANSITIONAL  FLOW 


Free  Stream 
Mach  No. 

Local 

Mach  No, 

Pnns/qe  (0  to  20  kHz) 

4 

3.7 

2.86  x  10~3  to  3.16  x  10-3 

8 

6.6 

1.85  x  10-3  to  2.70  X  10~3 

10 

8.1 

1.67  x  10-3  to  1.96  x  10*3 

Figure  30  represents  typical  normalized  spectra  (using  local  flow  properties)  for  transi¬ 
tional  flow.  Figure  51  gives  the  third  octave  pressure  distributions  at  various  points  in  the 
transition  region.  During  transition  onset  and  decay  it  was  noted  that  tin?  high  frequence 
portion  of  the  spectra  were  initially  affected.  At  the  point  of  peak  transition  levels  across 
the  complete  pressure  spectrum  (0  to  20  kHz)  are  increased  3  to  to  dL  when  compared  to 
tunnel  tare  noise  and  5  to  10  dll  when  compared  to  turbulent  flow. 

At  angles  of  attack  equal  to  tlie  half  cone  angle  of  tiie  test  vehicle,  separated  flow  occurs  on 
the  leeward  side.  Figure  32  represents  the  pressure  spectral  distribution  obtained  at  the 
aft  circumferential  array  at  a  model  attitude  of  *7.2*.  Sensor  15  represents  the  spectrum 
in  separated  flow  white  sensor  19  is  in  transitional  flow.  A  10  to  20  di>  variation  exists 
across  the  spectrum  which  could  result  in  significant  hoop  excitation,  due  to  Urn  harmonic 
distributions  of  pressure. 

Typical  base  pressure  spectra  at  Mach  4  are  shown  in  Figure  53,  measured  at  the  base 
center  for  various  angles  of  attack.  In  general  jauik  levels  were  measured  at  zero  angle  of 
attack.  As  Reynolds  number  increased  base  spectra  was  olxterved  to  Increase  propur- 
tionally.  Data  was  obtained  for  various  nose  radii,  however,  no  effects  were  observed  v«» 
base  sensors,  indicating  that  the  base  spectra  is  uir*fteeted  by  changes  in  bluno. » . 
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Fluctuating  Pressure  Distribution  Tunnel  Condition 
(Reoo  -  2.8  x  106,  M^-4,  Rn-0.0) 
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NORMALIZED  POWER  SPECTRAL  DENSITY 
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Figure  52.  Comparison  of  Pressure  Spectra  In  Circumferential  Direction 
(Separated  Flow  Leeward  Sensors),  M  =  8,  Re  =  2.8  x  I0G 
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DB  CRE  0-0002  HICROBfiR? 


CNE-TMIRO  OCTAVE  SAND  CENTER  FREQUENCY  IN  HZ 


Figure  53.  Comparison  of  Base  Flow  Pressure  Spectra  for  Various  v 
(Base  Center),  =  4,  Rew  =  2.2  x  10<> 
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CROSS  CORRELATION  FUNCTIONS 


*>. 


Space-time  correlation  functions  representing  turbulent,  transitional  and  separated  flow 
are  shown  in  Figure  54,  The  typical  decrease  and  broadening  of  peaks  for  increasing 
separation  distance  are  attributed  to  decay  of  high  frequency  turbulent  eddies  as  they 
progress  downstream.  The  only  significant  difference  between  turbulent  and  transitional 
flow  is  in  the  convective  velocity.  Correlation  functions  obtained  at  high  Mach  numbers 
(8,  10)  tend  to  be  broad  for  small  separation  distances.  This  effect  is  thought  to  be 
caused  by  tunnel  noise.  No  correlation  was  evident  for  base  flow. 

4.  CONVECTION  VELOCITY 


Convection  velocities  in  transitional  flow  were  approximately  0.6  times  that  for  turbulent  flow. 
Figure  55  compares  broad  band  convection  velocities  for  transitional  and  turbulent  flow 
for  various  separation  distances.  Turbulent  flow  convection  velocities  measured  are 
typical  to  that  of  other  experimenters  discussed  in  Section  4. 
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SECTION  VI 


PREDICTION  METHODS  REFINEMENT 


UtiJizing  existing  data  summarized  in  Section  II,  and  data  obtained  in  this  tunnel  program  as 
presented  in  Section  V,  Volumes  I  and  n,  Volume  n,  recommendations  for  the  aerodynamic 
properties  necessary  to  describe  the  aeroacoustic  loading  environment  are  given.  The 
criteria  for  recommended  methods  were  based  upon  the  most  reasonable  fit  of  experimental 
data  and  their  practicality  in  design  use. 

Attempts  to  develop  narrow-band  cross  correlation  functions  have  resulted  in  wide  scatter 
over  the  frequency  spectrum  of  interest.  This  apparent  erratic  behavior  in  the  narrow  band 
results  has  been  traced  to  limitations  in  the  Fast  Fourier  Transform  technique  used  to  ana¬ 
lyze  the  data.  While  not  totally  understood,  indications  are  that  the  limiting  process  utilized 
in  the  digital  program,  while  rigorously  correct,  is  not  compatible  with  the  finite  sample 
size  associated  with  the  data  of  this  study.  This  is  the  only  plausible  explanation  for  the 
apparent  nontypical  behavior  in  the  narrow-band  correlation  plots.  Hence,  this  information 
is  not  included  in  the  results  of  this  study. 

1.  ATTACHED  TURBULENT  BOUNDARY  LAYER  FLOW 


a.  Acoustic  Magnitude 

Figure  56  summarizes  equations  used  for  acoustic  magnitude,  as  proposed  by  various 
experimenters.  These  are  compared  with  test  data  including  data  obtained  in  the  present 
test  program.  Figure  57  presents  the  data  of  other  experimenters  including  present  data 
extrapolated  out  to  infinity,  for  purposes  of  comparison.  Lowson's  expression  compares 
well  with  available  data  below  Mach  5  and  tends  to  unde rp rediet  for  higher  Mach  numbers. 

This  is  due  to  the  available  test  data  trend  to  group  at  two  levels  (0,006  for  M  £  2  and  0.0015 
for  A  M  <  10).  In  order  to  be  representative  at  both  levels,  n  variable  Mach  number  co¬ 
efficient  in  the  denominator  of  Lowson's  expression  becomes  necessary,  Houbolt's  refined 
equation  accounts  for  this  by  use  of  a  coefficient,  re,  which  is  a  function  of  temperature  and 
local  velocity.  Table  VII  compares  the  present  tost  data  with  predictions  using  Houbolt's 
refined  expression  for  acoustic  magnitude.  Those  values  differ  from  those  presented  in 
Figure  56  in  that  actual  test  conditions  were  used  to  determine  the  recovery  factor  1*1  in  the 
expression  for  r0,  where  in  Figure  56  a  value  of  unity  for  the  recovery  factor  (ri)  was 
assumed  to  calculate  r0. 

Houbolt's  expression  (repeated  below)  for  acoustic  magnitude  is  therefore  recommended  for 
general  use. 
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Figure  56.  Comparison  of  Overall  Acoustic  Magnitude  Prediction  Methods  to  Experimental 


ce-nn  Square*  Pressure  Fluctuation  Vs  Mach  Number  for  Attached  Turbulent  Flow 
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Because  of  its  simplicity  and  close  fit  to  experimental  data  (M  <  5)  Lowson's  equation 
(Equation  1)  can  also  be  applied.  Above  Mach  5,  the  equation  should  be  adjusted  to: 
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TABLE  VII,  PREDICTION  VS.  TEST  ACOUSTIC  MAGNITUDE  (TURBULENT  FLOW) 


Free  St:  earn 

Maoh  No. 

P  /q  x  103 
rms  M  e 

Test  Data* 

Predicted  (Equation  3) 

4 

1.53  to  1,9 

2.5 

8 

0.835  to  1.58 

1*5 

10 

1.08  to  2.04 

1.2 

♦Tost  data  adjusted  to  overall  value  for  comparison  with  Equation  11  and  previous  data, 
b.  Acoustic  Power  Spectral  Density 

The  criteria  used  in  defining  an  expression  for  power  spectral  density  function  wore 
that  the  integral  of  the  function  is  equal  to  the  root  mean  pressure  squared  (  f*  $  (f)  df  « 
p^rms),  the  function  remain  finite,  and  that  the  function  be  symmetric  about** ero,  a  require¬ 
ment  due  to  its  Fourier  transform  relationship  with  the  autocorrelation  function.  Figures  58, 
59,  and  60  compare  Houbolt'a  and  Lowson's  equations.  Both  expressions  compare  well  with 
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NORMALIZED  POWER  SPECTRAL  DENSITY 


U» 


Figure  58.  Normalized  Power  Spectral  Density  Distribution  for  Turbulent  Flow  M  w 


TUNNEL  TARE  ELIMINATED 


Figure  59.  Normalized  Power  Sjpoctral  Density  Distribution  for  Turbulent  Flow  M 


NORMALIZED  POWER  SPECTRAL  DENSITY 


Figure  GO.  Normalized  Power  Spectral  Density  Distribution  for  Turbulent  How  M 


present  data  (except  at  Mach  10)  with  Houbolt's  (Equation  10)  expression  being  reasonably 
conservative  in  the  low  frequency  range  of  primary  interest.  Since  Lowson's  spectral  equa¬ 
tion  (Equation  11)  is  defined  in  terms  of  boundary  layer  thickness  ($)  it  was  necessary  to 
express  6  in  terms  of  the  boundary  layer  displacement  thickness  given  in  Table  III  of  Section 
V .  Present  data  (M  oo  =  4, 8)  are  also  shown  in  Figure  61  with  data  obtained  by  other  experi¬ 
menters  .  Good  correlation  was  obtained . 

Of  the  relationships  discussed  in  Section  IV,  Houbolt’s  expression  (Equation  10) 
given  below  is  recommended  for  general  use  because  of  its  ease  of  application  and  degree 
of  conservatism. 
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c.  Cross-Correlation  Functions 

Due  to  the  inability  to  derive  representative  narrow-band  cross  correlation  func¬ 
tion^  from  the  wind  tunnel  data,  recourse  was  made  to  the  literature  for  typical  expressions 
icable  to  supersonic  flow.  Available  data  in  the  Mach  0.3  to  2.5  region  indicate  the 
form  of  tim  functions  to  be  as  given  in  the  following  equations: 
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Lowson*“*  and  Kobertson*1*  using  data  obtained  by  Bull*4*  and  Chyu  and  Manly*9* 
defined  the  coefficients  for  the  above  equations  to  be: 


A„(»),  u>)  *  exp  (-0.72J  a  |~p)  exp  (-2.Q  jnj) 

U«  & 


(lb) 


h  (*,w)=exp  (-O.ljl  j-- )  exp  (-0.27  hi) 

c  5 


(19) 
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The  above  expressions  (Equations  18  and  19)  may  also  be  written  in  terms  of  the  boundary 
layer  displacement  thickness  parameter  as: 


At,(i?,w)  =  exp  (-0.72|fi|^-)  f  0.3  +  0.7  exp  (  -O.jlhl  )] 

UC  L  5  *  J 

A*  (£,co)  =  exp  (-0 . 1  Jli w )  exp  (-0.034  lit  ) 


(20) 


(21) 


Equations  16  through  21  are  recommended  in  defining  the  cross-power  spectra  for  turbulent 
flow. 

d.  Convection  Velocity 

Broad  band  convection  velocities  obtained  in  the  present  test  program  (M  »  =  4  to 
10)  were  typical  of  those  observed  by  other  experimenters  (M*.  =  0.3  to  2.5),  therefore 
Equation  24  developed  by  Lowson^  for  the  narrow  band  convection  velocity  in  turbulent 
flow  is  considered  representative. 

U 

=  0.075  +0.3  exp  (-O.U^r  )  -°*25  H*2  £  )  (24) 

UoO  r~ 


2.  TRANSITIONAL  BOUNDARY  LAYER  FLOW 

Figure  46  gave  typical  acoustic  definition  of  the  fully  developed  turbulent  boundary  layer 
(TBL)  and  transitional  aerodynamic  flow  regions  for  the  various  wind  tunnel  conditions  and 
model  orientations  investigated  in  this  program.  Curves  which  connect  the  data  points  pre¬ 
sented  therein  were  generated  in  conjunction  with  aerodynamic  prediction  of  the  transition 
region  and  represent  the  estimated  shape  of  the  sound  pressure  level  (SPL)  profile.  For 
the  purpose  of  evaluating  acoustic  data  ucar  the  peak  of  transition,  analysis  was  restricted 
to  those  specific  cases  where  the  transducer  registered  within  2  d!3  of  the  peak  transitional 
fluctuating  pressure  level .  Once  normalized  values  of  the  acoustic  magnitude  (Prms)  and 
power  spectral  density  functions  ($(f) )  had  been  obtained  for  these  cases,  the  multiplicative 
factors  given  in  Table  VIII  were  utilized  to  calculate  values  representative  of  peak  transition. 
When  the  data  point  location  on  the  SPL  profile  is  at  the  estimated  peak  no  adjustment  was 
required . 


Note:  Contributions  from  tunnel  tare  noise  (Sensor  2)  are  sufficiently  lower  than  measured 
values  of  transitional  spectra  and  acoustic  magnitude  that  these  noise  contributions  are  con¬ 
sidered  negligible  in  the  analysis  of  transitional  data  which  follows. 
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a.  Acoustic  Intensity 


The  value  of  Prms  (0  to  20  kHz)  for  a  given  condition  is  typically  presented  by 
normalizing  with  respect  to  either  free  stream  or  local  dynamic  pressure  (qM  ,  qe)  It 
was  determined  that  a  tighter  collapse  of  the  present  data  resulted  from  a  qe  normalization. 
However,  for  comparison  with  other  transition  data,  the  measured  values  depicted  in 
Figure  62  were  normalized  (after  adjustment  per  Table  VIII)  with  respect  to  the  free  stream 
dynamic  pressure.  Using  qe  as  the  normalizing  parameter,  Figure  63  reveals  how  the 
measured  values  of  Prmg/qe  near  peak  transition  exceed  those  determined  for  turbulent 
boundary  layer  flow.  In  order  to  develop  a  refined  acoustic  intensity  prediction  technique 
for  transitional  pressure  data,  an  attempt  was  made  to  extend  the  frequency  range  of  anal¬ 
ysis  out  to  infinity.  This  was  accomplished  by  employing  the  equation  for  acoustic  power 
spectral  density  derived  in  Section  VI. 2. b  to  fit  the  spectral  data,  i.e. 

xo  - — p - ; - TT  <51* 

i  +|c(mq)  6  f/ue  j 

where  C(Mq)  is  given  by  Equation  50. 

Therefore  integrating  both  sides  of  Equation  51  from  zero  to  Infinity  yields  the  following 
formula  for  the  mean  square  aooustio  pressure  in  terms  of  tho  local  dynamic  pressure  qy: 

2  3.14  x  Iff4  qft2 

RMS  =  CfM~)  (45) 

Tho  plotted  results  are  shown  in  Figure  64  which  indicates  versus  local  Mach  number 

as  predicted  from  the  spectral  data  using  Equations  51  and  43. 

Values  for  transitional  acoustic  magnitude  given  hi  tills  figure  reveal  a  generally 
shallow  data  trend  for  4.0  <  <  10 .0 .  The  equation  for  rma  pressure  which  describes 

those  plotted  values  in  the  Mach  number  range  investigated  is 


(3.7  <,  Me  £  8.1) 


(46) 


The  above  expression  is  quite  similar  in  form  to  Lowson's  equation  for  intensity  of  fully 
turbulent  flow  -  0,006/(1  +  0, 14  M2C).  An  assumption  that  these  two  equations  are 

both  applicable  at  low  Mach  numbers  would  imply  a  cross-over  of  Ppf,jtj/qc  values  fer 
transitional  and  fully  turbulent  aerodynamic  conditions  as  Me'“*‘0.  However,  it  is  evident 
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from  Figure  56  that  normalized  values  of  fully  turbulent  acoustic  magnitude  -follow  different 
trends  for  Mach  numbers  less  than  and  greater  than  approximately  four.  Therefore,  ex¬ 
pressions  similar  to  Lowson's  must  be  qualified  with  the  applicable  Mach  number  range. 

This  suggests  that  the  same  trend  holds  for  transitional  data,  such  that  an  equation  some¬ 
what  different  from  Equation  46  would  have  to  be  employed  for  transitional  intensity  at  low 
Mach  number.  Quantitative  effects  of  transition  at  subsonic  and  transonic  velocities  should 
therefore  be  the  subject  of  further  investigation.  However,  for  high-beta  re-entry  vehicles 
transition  is  associated  with  high  Mach  numbers  and  expressions  derived  in  this  section 
defining  the  acoustic  properties  of  transitional  flow  are  directly  applicable  to  vehicles  of 
this  type . 

b .  Acoustic  Power  Spectral  Density 

For  those  cases  of  transitional  flow  listed  in  Table  VIII,  the  one -third  octave  band 
SPLdata  were  employed  to  calculate  spectral  values  versus  frequency  out  to  20,000  Hz. 
These  values  of  $  were  normalized  to  the  local  aerodynamic  TBL  parameters  Ue,  *s  *,  and 
qe.  Spectral  values  are  shown  plotted  in  Figaros  65  to  67  as  functions  of  the  non-dimensional 
frequency  parameter  P  «  f  5  */Ue  (Strouhnl  number)  for  the  various  test  Mach  numbers. 

No  clearly  identifiable  spectral  shape  difference  was  noted  between  apparent  peak  transition 
PSD's  and  those  estimated  to  be  within  2  dB  of  the  apparent  peak.  Therefore  levels  of  <£> 
were  adjusted  uniformly  throughout  the  20  kHz  frequency  range  using  the  factors  given  in 
the  preceding  section.  From  Figures  65,  66  and  67,  it  can  be  seen  that  a  reasonable 
collapse  of  the  resulting  normalized  data  is  obtained  and  Unit  a  definite  trend  in  the  data  Is 
apparent.  It  is  postulated  that  the  data  curl  up  at  lUgh  frequency  exhibited  by  Mach -  -l 
spectra  results  from  an  overcorreetion  for  sensor  size  effects  in  that  limited  frequency 
range. 


A  comparison  of  those  figures  reveals  that  normalized  spectral  values  for  a  Strouhul 
number  less  than  10“2  tend  to  be  independent  of  Mach  number,  but  that  the  nondi  mens  tonal 
cut  off  frequency  f  '  is  less  for  M  ^  •  10  than  for  M«*  6  and  l»  These  facts  are  taken  into 

account  below*  in  developing  prediction  equations  for  Ph&j$Ai0  sad  . 

Figure  68  is  a  plot  of  a  typical  4  PSO  from  tl»c  present  investigation  super¬ 
imposed  on  the  spectral  data  from  Figure  15.  As  mentioned  in  Section  111,  the  high  levels 
seen  at  low  frequency  in  tiic  previously  obtained  data  may  be  due  to  tunnel  noise  and  elec¬ 
trical  effects . 

Previous  techniques  for  determining  predicted  levels  for  a  transition  environment 
involved  using  the  turbulent  boundary*  layer  acoustic  prediction  corresponding  to  given 
aerodynamic  flow  conditions  as  a  baseline.  Then  die  limited  transition  data  available  was 
employed  to  scale  up  the  TBL  spectral  levels  to  estimated  peak  transition  values.  There¬ 
fore  it  is  important  to  use  a  comparison  of  TBL  and  transitional  spectra  from  the  present 
data  to  verify  and/or  refine  this  prediction  procedure.  As  an  example.  Figure  69  depicts 
normalized  transitional  PSD’s  dkmsars  11,  13  adjusted  to  apparent  peak)  and  compares 
these  data  to  turbulent  boundary  layer  PSD’s  for  Hun  20  at  M  4.  Normalized  spectra 
for  this  ease  exceed  TBL  values  by  approximately  1  decade  (10  dB) . 


n  mo 


iy 

ZT 


Figure  66 .  Normalized  Power  Spectral  Density  Distribution 
for  Transitional  Flow  M*. «  8 
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It  is  the  intent  of  this  section  to  derive  a  method  for  prediction  of  transitional 
quantities  more  precise  than  the  previous  indirect  approach,  but  using  a  formal  development 
quite  similar  to  that  of  turbulent  acoustic  analysis.  Attempts  at  generating  an  empirical 
equation  to  describe  TBL  spectra  have  generally  assumed  that.the  PSD  level  as  a  function  of 
frequency  has  a  form  similar  to 


*(f) 


K,  P* 

1  RMS 


i  +  <k2  f) 


(47) 


2 

such  that  <J>(f)— for  low  frequencies.  Values  of  are  dependent  on  aerodynamic 

parameters  of  velocity,  boundary  layer  displacement  thickness  and  dynamic  pressure . 
Examination  of  Flgur  s  65  and  67  suggests  that  peak  transition  spectral  data  may  be  described 
by  an  equation  of  similar  form,  with  p|!  replaced  by  qQ2.  Therefore  defining  nondimen- 
sional  quantities 


=  4>U  A  * 
e 


and  f  =  f5  */U 

e 


we  get 


*'(£') 


(46) 


where  the  character  of  the  undetermined  coefficients  are  defined  by  a  best  fit  of  the  data. 

An  equation  for  which  provides  a  relatively  good  fit  to  the  normalized  data  is 
the  following} 


$'(f’)  “ 


2.00  x  IQ**4 

1  +  [c(Me)f]  2 


(49) 


where 

2  2 

C(M  )-i8.6  (  1  +0.013  hi  ) 

0  6 

Figures  65,  66  and  67  present  Equation  49  superimposed  on  the  normalized  spectra. 
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(50) 


In  terms  of  absolute  quantities.  Equation  49  becomes 


9(f)  = 


2.0  x  10~4 5  *  q  2 /U. 

G  v 


l  + 


*  ”1* 

C(M  )8  f/U  ‘ 
6  6 


(51) 


The  particular  form  of  C(Me)  in  Equations  49  and  51  was  chosen  for  compatibility  with 
Equation  46,  such  that  performing  an  integration  of  Equation  51  from  0  to  infinity  and  taking 
the  square  root  yields 


P  /q  = 
RMS/  e 


1.78  x  10 


vw 


-2 


0.0041 


1  +  0.013  M 

e 


which  is  identically  Equation  46 . 

Further  insight  into  the  implications  of  Equation  51  can  be  gained  by  examining 
the  identity 

0  ,*>  2.0xl0"4q2  w/2 

P«MS  \i  * - - 


C<M0) 


(521 


which  Implies 


2. Ox  10-4  q2 
e 


2P&m£,mo> 


(5.S) 


Substituting  this  into  Equation  51,  the  expression  for  transitional  acoustic  spectra  in  terms 
of  TBL  aerodynamic  parameters  becomes 


C  (MJ 


*(D 


2  IT 


2 


1  +  2* 


C(Mc) 

2* 


(54) 


6  f/U 


Now  define  a  fictitious  transitional  flow  length  parameter  6  (analogous  to  the  turbulent 
boundary  layer  displacement  thickness),  by  1 


C  (Me)  * 
2»  6 


3.7  ^Me  5  8.1 


(35) 
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Equation  54  hence  becomes 


4  <««  /U  )  PZ 

«f)=  - F - Ti  3.7  <  Me<6.1  (56) 

l+|^f5f/Ue  J 

which  is  Equation  11  for  TBL  spectra  expressed  in  the  frequency  domain  except  6  *  and  Uc 
have  been  respectively  replaced  with  Sj.  and  U  .  The  similarity  in  form  further  ensures 
that  both  sides  of  Equation  56  integrate  identically  to  P^g*  Therefore  it  is  postulated  herein 
that  transitional  spectra  in  the  range  3. 7  <  Me  8. 1  can  be  defined  using  a  prediction  tech¬ 
nique  completely  analogous  to  fully  turbulent  procedures,  using  Equations  55  and  56.  Figure 
70  describes  the  dependence  of  the  ratio  6^  /<5  *  on  local  Mach  number. 

c.  Cross  Correlation  Functions 

As  previously  discussed,  narrow  band  cross  correlation  functions  were  not  obtained. 
Since  the  spectral  distribution  for  transitional  flow  measured  in  the  present  test  program  was 
similar  to  that  of  turbulent  flow  it  appears  reasonable  that  the  Cross  Correlation  Function 
remain  the  same.  Until  further  resolution  of  this  area  is  made,  it  is  recommended  that  the 


Figure  70.  Functional  Relationship  of  Transitional  and  Fully  Turbuient 
Characteristic  Length  Parameters 
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correlation  functions  for  fully  developed  turbulent  flow  be  used  as  a  first  approximation. 
Equations  20  and  21  are  therefore  recommended  to  define  the  cross  correlation  coefficients 
for  transitional  flow. 


(20) 


(21) 


d.  Convection  Velocity 

Broad  band  convection  velocities  for  transitional  flow  determined  from  data  obtained 
in  the  present  test  program  were  approximately  0.6  of  that  in  turbulent  flow;  therefore,  the 
following  expression  is  recommended  for  transitional  flow 


U 

— ~  =  0.045  +  0.18  exp  ( 


-0.11 


)  -  0. 15  exp  ( 


(57) 


3.  BASE  FLOW 


a.  Acoustic  Magnitude 

Houbolt^0*  was  the  only  source  found  which  defined  expressions  for  the  acoustic 
environment  acting  on  the  base  of  a  ro-cntx'y  vehicle.  Based  upon  engineering  logic  and 
the  knowledge  of  turbulent  flow,  the  following  expression  was  domed  relating  base  pressure 
fluctuations  to  the  static  pressure 


P 

rms. 

b 


0.01  MjJ 

1+ 


P, 


b 


(38) 


(18) 

This  relationship  agroed  well  with  data  obtained  by  Eld  red'  1  for  subsonic  flow 
but  as  stated  in  Houbolt’a  report,  verification  of  the  coefficients  in  the  above  equation  was 
in  order  for  supersonic  flow.  Figure  71  compares  Equation  38  to  data  obtained  in  this 
study  at  zero  angle  of  attack  for  Mach  numbers  associated  with  the  wake  boundary  of  6.2 
and  13,6  (M«,=  4  and  10  respectively).  Predicted  levels  are  lower  than  measured  values 
indicating  that  adjustments  to  the  denominator  In  the  expression  given  by  equation  38  are 
required.  No  adjustment  to  the  numerator  coefficient  is  required  since  it  agrees  with 


Figure  71.  Overall  Acoustic  Intensity  for  Base  Flow  (Prediction  vs  Test) 


Eldrod’s  subsonic  results.  Using  a  coefficient  of  0,04  in  the  denominator  as  given  in 
Equation  58  agrees  well  with  data.  Figure  71  is  a  plot  of  Equation  68  and  the  tunnel  data. 


P 

rmsb 


0.01 

1  +  0.04 


(58) 


At  large  Mach  numbers,  Equation  68  suggests  that  the  rms  pressure  level  is  dependent 
simply  on  the  base  static  pressure,  reaching  an  asymptotic  value  of  0,25  P^. 

b.  Acoustic  Powor  Spectral  Density  and  Correlation  Function 

In  order  to  arrive  at  the  power  spectral  density  for  base  flow  Houbolt  assumed  a 
correlation  function  of  the  form 


R  (r)  o  P; 


-  (VT/L)t 


rasb 


C39) 


where  VT  was  an  effective  transport  velocity  and  L  the  scale  of  turbulence.  Assuming  that 
the  transport  velocity  was  related  to  the  transverse  component  of  the  wake  boundary  velocity 
(VT  =  Vb  sin  0  where  6  is  the  half  cone  angle  of  the  wake  cone)  and  the  scale  of  turbulence 
(L)  proportional  to  the  base  radius,  Equation  59  was  developed  for  the  base  power  spectral 
density: 


<t>  (u>)  = 


2  P  L 
rms, 
b 


n  V„ 


2  P 

rms,  r, 
_ b  b 

jtU,  sine 
b 


/_VLV 

K  sin7 


(59) 


1  + 


Equation  59  is  compared  with  normalized  test  data  (Meo=  4,  Rq  -  1.4  x  106, 

2.2  x  10®,  2. 8  x  10®)  in  Figure  72,  It  is  evident  that  this  equation  is  conservative  at  low 
Strouhal  numbers.  Figure  73  compares  the  same  Mach  4  data  with  Equation  60  which  assumes 
that  the  transport  velocity  is  equal  to  the  wake  boundary  velocity.  Good  agreement  is 
attained  with  test  data  therefore  Equation  60  and  its  associated  correlation  function.  Equation 
61  are  recommended  to  define  acoustic  environments  for  base  flow.  Figure  74  compare  the 
recommended  equation  to  data  obtained  at  =  10  reasonably  good  correlation  is  attained: 
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APPENDS  X 

ENGINEERING  RELATIONS  /GRAPHS  FOR  THE  DETERMINATION 
OF  LOCAL  FLOW  PROPERTIES  ON  SHARP  BODIES 

I 

A.  MARTELLUCCI 


Within  the  body  of  this  report,  the  measured  acoustic  data  have  been  normalized  with 
various  aerodynamic  flow  properties  such  as  the  boundary  layer  edge  dynamic  pressure, 
s  1/2  (pe  the  boundary  layer  displacement  thickness,  6*,  the  local  static 
pressure  etc.  It  is  the  purpose  of  this  appendix  to  provide  the  reader,  who 

may  not  have  ready  access  to  viscous  flow  computer  programs  or  who  may  be  unfamiliar 
with  the  equations/ techniques  necessary  to  compute  the  viscous  properties,  with  sim¬ 
plified  equations  and  charts  with  which  to  determine  some  of  the  more  fundamental 
local  flow  properties. 


In  general,  for  blunted  bodies,  the  determination  of  the  properties  within  the  viscous 
layer  involves  ;he  numerical  solution  of  the  boundary  layer  equations  using,  for  example 
a  finite  difference  scheme.  Another  common  approach  involves  casting  the  equations 
in  an  integral  form  resulting  in  integro-differential  equations  which  are  solved  by  Inte¬ 
gration  along  the  body  surface  to  the  station(s)  in  question.  As  a  result,  computer  pro¬ 
grams  have  been  formulated  to  solve  these  systems  of  equations.  Programs  exist  for 
the  solution  of  the  laminar  and  turbulent  boundary  layer  flows,  however,  when  the 
boundary  layer  on  the  body  is  laminar  -  transitional-turbulent,  special  devices  or  tech¬ 
niques  must  be  employed  to  solve  this  class  of  problem.  In  the  finite  difference  solution 
scheme,  in  which  a  forward  marching  solution  is  employed,  one  customarily  lets  the 
viscosity  vary  in  the  transitional  boundary  layer  region  (which  must  be  known  a  priori) 
from  the  molecular  value  at  transition  onset  to  the  turbulent  effective  value  at  tiro  end 
of  transition  using,  say,  a  gaussian  distribution. 

The  integral  form  of  the  viscous  layer  equations  are  more  commonly  used  in  vehicle 
design.  In  this  type  of  scheme,  it  is  quite  common  to  treat  the  transitional  region  as 
a  point,  that  is,  the  viscous  layer  is  treated  as  laminar  to  this  point  ami  downstream 
of  this  point  the  turbulent  flow  relations  aro  used.  To  pass  through  this  region  one 
must  establish  the  constants  of  integration  for  the  turbulent  solution  (i.e. ,  the  so-called 
"effective  origin"  concept).  This  is  customarily  performed  by  equating  the  laminar 
momentum  thickness,  0^  ,  'o  turbulent  valuo,  0 ^  ,  thereby  establishing  an  artificial 
origin  to  the  turbulent  solution  such  that  at  that  particular  station  of  matching  (usually 
transition  onset)  tire  value  of  °  l  »  6  -p  .  The  remaining  turbulent  properties  can  then 
be  determined.  This  information  is  by  way  of  background.  I  will  attempt  to  present 
simplified  laminar  and  turbulent  boundary  layer  equations  valid  for  pointed  cones  or 
sharp  flat  plates.  Furthermore,  graphs  will  be  presented  from  which  the  inviscid  flow 
(boundary  layer  edge)  properties  oan  bo  determined.  It  should  again  bo  stressed  that 
the  following  information  is  valid  only  for  sharp  leading  edge  bodies.  Strictly  speaking, 
for  blunted  bodios  the  integral  equations  must  be  solved.  Thus  tiro  following  information 
provides  a  means  of  estimpMng  approximate  values  or  estimates  of  parameters  which 
are  valid  for  blunt  body  stations  far  from  the  leading  edge. 
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INVISCID  FLOW 


In  this  section,  a  series  of  graphs  will  be  presented  which  can  be  used  to  define  the 
inviscij  perfect  gas  (v  =  1.40)  flow  properties  on  pointed  cones  and  wedges.  Listed 
below  are  the  Figure  numbers  and  the  flow  properties  presented  in  each  chart. 


PROPERTY 

Wall  Pressure  Ratio,  Pw/P  =  P J  P  ,  vs  M  for 
various  9C 

Edge  Temperature  Ratio,  Tt/T  vs  M  for  various 

Ac 

Edge  Mach  number,  M„  vs  M  for  various 

C  00  ^ 

Local  unit  Reynolds  number  ratio,  Ry/R  '  s 
for  various  M 


With  the  free  stream  conditions  known,  the  local  dynamic  pressure  on  the  body,  qy, 
can  be  determined  from  the  follow  big: 


io  a  1/2  o0  =  x  (-^r) 


where  ?/2  =  0. 7 

The  relation  for  the  local  velocity  is: 


u„  -  M„ », 


*  M  «  y  11T  •  45. 01  M 


The  local  wetted  length  Reynolds  number  at  a  distance,  s,  along  the  body  surface  can 
be  determined  from  Figure  A-4  or  A- 9  and  the  following; 


(v\) 

unit  Reynolds  number, 


R„  S  where  11„  gP  u  /  u  is  the  free  stream 

®  ®  te 


BOUNDARY  LAYER  EQUATIONS  -  Compressible  Flow/Ideal  Gas 

Viscous  layer  equations  will  be  presented  here  to  enable  the  reader  to  calculate  the 
momentum  or  displacement  thickness  for  a  cone  or  wedge  with  either  all  laminar,  all 
turbulent,  or  a  laminar-turbulent  boundary  layer  (using  the  matched  momentum  thick¬ 
ness  approach). 

LAMINAR  FLOW 


Boundary  Layer  velocity  thickness : 


K  (T*A,) 


where 


K  =  0  flat  plate,  K  =  1  cone 
Reference  Temperature: 

T*  =  0. 5  T  +  (0. 5  -  0,22  r)  T  +  0. 22  r  Tn 

W  U  U 

a 

Recovory  Temperature : 

Tr  =  r  T0  +  (1-r)  T0 


Viscosity: 


3/2 


u  “  2.270 


T 


x  10' 


•8 


LB-Soc 


T+ 108,0  ft* 

r  =  recovery  factor  »  0. 848  (laminar) 
recovory  factor  »  (.  896)  (turbulent) 
and  Tu,  =  model  wall  temperature 
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(4) 


(5) 


(0) 


(7) 


The  ratio  of  the  boundary  layer  displacement  to  velocity  thickness  is  given  by : 


v« 


j] 

S?  jt- 1 


6L*  ue2/2he  +  3,36  T\/Te  _0*376 


6L  ue2/2he  +  3.36TyTe  +  5.79 
where 


(8) 


he  "  Cp  Te 


(9) 


ft2 


Cp  =  6006 


sec2  °R 


(for  air) 


The  ratio  of  the  boundary  layer  momentum  to  velocity  thickness  is  given  by: 

_  ,  0. 18 

fiL  /  T* 

—  =  0,125 


(10) 


Through  the  use  of  equations  (4)  through  (10),  tho  local  viscous  parameters  f  . ,  * 

and  0  *  can  be  determined,  ^  L 

TURBULENT  FLOW 


Ihe  momentum  thickness  can  be  expressed  as: 

0.037  c  S 

®?  °  (11) 
l.i)UK  Ho  °-2 

o 

where 


Ihe  boundary  layer  displacement  thickness  can  be  obtained  from  the  follow  big  relations, 
in  conjunction  with  equation  (11): 


A-4 


Thus,  through  the  use  of  Figures  (A-l)  through  (A- 9)  in  conjunction  with  equations  (1) 
through  (13)  those  local  properties  employed  in  the  normalization  of  the  forebody 
acoustic  measurements  can  be  estimated. 
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BASE  PRESSURE  CORRELATION 


The  RMS  acoustic  pressure  on  the  base  of  conical  bodies  is  normalized  with  the  base 
pressure.  When  the  flow  at  the  cone  base  is  transitional- turbulent,  experimental  data 
have  shown  that  the  pressure  across  the  base  is  approximately  constant.  A  correla¬ 
tion  of  available  flight  and  wind  tunnel  cone  data,  results  in  the  line  shown  in  Figure 
A-10,  where  PjJ/Pw  is  plotted  against  Me»  P  and  M  can  be  determined  from  Figures 
A-l  and  A-3,  respectively. 

FLIGHT  EXAMPLE  ON  THE  USE  OF  THE  APPENDIX: 

Let  us  assume  that  one  wishes  to  determine  the  acoustic  magnitude  levels  and  the 
power  spectra  for  a  conical  re-entry  vehicle  at  an  altitude  of  10,000  feet.  The  vehicle 
in  this  case  is  6°  half  angle  cone  which  is  6  feet  long  and  has  a  nose  to  base  bluntness 
ratio  R»j/Rg  <  .02  (hence  a  sharp  cone  approximation  for  this  case  will  be  valid). 

Let  us  farther  assume  that  at  H  =  10  Kft,  Mw  =  15  and  Tw  =  3000°R. 

Fi'om  standard  altitude  tables  (e.  g.  1962)  one  can  establish 

Pm  =  1455. 6  PSFa 

To,  =  483°R 

Tq  =  22,555°R* 

CO 

RVft  =  8.025  xl07 

Thus  the  goomotry  and  free  stream  conditions  aro  set,  and  aro  summarized  below: 


♦It  should  be  noted  here  that  for  the  flight  case,  the  freestreum  stagnation  properties 
must  be  obtained  from  a  Mollior  Diagram.  For  simplicity  here,  ideal  gas  relations 
were  used. 


The  local  inviscid  flow  properties  can  be  determined  from  Figures  A-l  to  A- 4  and  are 
listed  below: 


Fig.  A-l 

P  /P  =5.0 
w  » 

Pw=  7278  psf 

Fig.  A-2 

T  /T  =  1.73 

T  =  835. 6°R 

6  co 

e 

Fig.  A-3 

M  =11.4 
e 

Fig.  A- 4 

R  /R  =  1.86 
e  e 

R  /ft  =  1.49  x  10 
e 

» 


q  =  y/2  M  2  P  =  662,000  psf 
e  e  w 

Using  the  acoustic  intensity  prediction  suggested  in  the  body  of  this  report  for  fully 
developed  turbulent  flow,  i,  e. , 

P  /q  =  0. 002  / 1  1  +  0. 02  M  2  1  =  0. 000555 
rms  ne  /  \  e  / 

which  corresponds  to  179  dB. 

For  the  fully  developed  turbulent  flow  case  presented  here,  one  can  readily  compute 
the  displacement  thickness  using  the  equations  presented.  The  procedure  would  be  as 
follows: 


Given  T  ,  T  .  and  the  freestream  conditions  one  can  determine  T*  and  T  from 
equations  (5)  and  (6),  respectively.  That  is 

T*  =  0. 5(3000)  +  5  -  . 22(.  896)J  835.  6  +  0. 22(.  896)  (22,555) 

T*  =  6,199°R 


T  =  0. 896  (22 , 555)  +  (1  -  .  896)  835, 6 

T  =  20,296°R 
r  * 


A-7 


•,  9  *  ’  ^  r  -  y  _  +?[  * -  rrgg  1  yrwy"  :  *  '■***  •»  jr  •.  - :  -  -  »xr  —  ■-  •  _-«* | «_ 


Before  computing  (L  from  equation  (11),  one  must  first  determine 

/,,*\  °-2  /m-\  0.8 


-  ©  (*) 

r  /16833  \  3/2  7835.6  +  198.  6\~]  0.2  l~835.6~]  0. 
y  835^6  y  ^6199  +  198.6  )  6199 


€  =  0.255 

Thus  from  equation  (11)  at  S  =  6' 
„  0.037  (.255)  (6) 


T  1.913  [1.49? 

=  0. 000479  Ft. 
Then  from  equation  (12) 


49xl08)  (6)10,  2 


=  0.  000479 


-i+^| 


«■«** 


i*l] 


6*t  =  0.  0309  Ft. 


This  value  would  be  used  in  the  acoustic  power  spectral  density  relation  /f  (w)  as 
defined  in  the  body  of  the  report. 

This  numerical  example  serves  as  an  illustration  in  determining  the  normalization 
parameters  utilized  in  the  body  of  the  report.  It  should  be  recognized  that  for  the 
flight  caso  real  gas  flow  tables  should  be  employed.  For  wind  tunnel  cases,  ideal  gas 
flow  relations  and  tables  are  generally  adequate. 


L 
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Figure  A- 


Figure  A-2 


A- 10 


Figure  A-3 
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